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Spin relaxation due to the D'yakonov-Perel' mechanism is intimately related with the spin split- 
ting of the electronic states. We determine the spin relaxation rates from anisotropic spin split- 
tings of electron subbands in n-(OOl) zinc-blende semiconductor quantum structures calculated self- 
consistently in the multi-band envelope function approach. The giant anisotropy of spin relaxation 
rates found for different spin-components in the (001) plane can be ascribed to the interplay between 
the bulk and quantum well inversion asymmetry. One of the in-plane relaxation rates may exhibit 
a striking nonmonotonous dependence on the carrier density. 
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I. INTRODUCTION 

The spin degree of freedom of electrons in solids has 
recently attracted much interesfeiiSiiiii in the perspective 
of spintronic devices, which require long spin relaxation 
times. Thus, a quantitative understanding of the depen- 
dence of spin-relaxation times on the system parameters 
is needed for the desired engineering of systems with the 
required properties as well as for the interpretation of 
measured values iSiS*^ It is generally accepted^ that in 
bulk zinc-blende type semiconductors and in quantum 
well (QW) structures based on these materials electronic 
spin relaxation is governed by the D'yakonov-Perel'iiiii^ 
(DP) and the Bir-Aronov-Pikus^- mechanisms. The for- 
mer mechanism becomes dominant in n-doped (001)- 
grown QW systems, which we consider here. The DP 
mechanism is intimately related with the spin splitting 
of the electronic subbands caused by the lacking inver- 
sion symmetry. "'^'^ It can result from the bulk crystal struc- 
ture (bulk inversion asymmetry, BIAi^) but also from the 
QW structure (structural inversion asymmetry, SIA^^). 
Making use of this relation, Averkiev et dlMM. employ 
a simple variational calculation for a triangular poten- 
tial model of n-doped GaAs heterostructures and include 
BIA and SIA terms as a perturbation to calculate spin- 
splitting and DP spin-relaxation rates. The simplicity 
of this approach raises the question of the reliability of 
these results in view of a more rigorous calculation, as 
will be presented in this paper. 

In this contribution, we determine the electron spin 
relaxation rates for longitudinal (in growth direction) 
and transverse (in-plane) spin components in asymmetric 
n-doped (OOl)-grown zinc-blende semiconductor quan- 
tum structures using the anisotropic spin splitting ob- 
tained from self-consistent calculations in the multi-band 
envelope- function approximation (EFA).^^ We investi- 
gate the dependence of these rates on the QW width 
and carrier concentration for the material systems Al- 



GaAs/GaAs and AlGaSb/InAs, and compare our results 
for the former system with those of Ref. [Ij. 



II. DP MECHANISM AND SPIN SPLITTING 

The DP spin relaxation is described in a single band 
approach by considering the time evolution of the elec- 
tron spin-density vector 



S = (tr(pfccr)) = ^tr(pfecr). 



(1) 



where k is the electron wave vector, <j the vector of 
Pauli spin matrices, and pk the 2x2 electron spin- 
density matrix. The dynamics of pk is ruled by a 
Bloch equation containing the spin-orbit coupling Hso = 
her ■ J1(A;) and the (spin independent) momentum scat- 
tering, e.g., with phonons, (nonmagnetic) impurities or 
other electronsiSiSS characterized by a momentum scat- 
tering time Tp. The spin-orbit coupling is analogous to a 
Zeeman term, but with an effective magnetic field r2(fc) 
that depends on the underlying material and geometry 
and on the electron wave vector k. The vector tr(pfc<T) 
precesses about $!(&) which eventually results in spin re- 
laxation. However, momentum scattering changes k and 
thus the direction and magnitude of Vl(k) felt by the elec- 
trons. Therefore, frequent scattering (on the time scale 
of |J1(A;)|~^) suppresses the precession and consequently 
the spin relaxation. This is the motional-narrowing be- 
havior that is typical for the DP mechanism, according 
to which the spin relaxation rate t^T^ (x 

Due to time reversal invariance we have $7(0) = 0. At 
finite k the spin-orbit coupling Hso causes a spin splitting 
2h\n{k)\ of the subband structure. It exists already in 
bulk semiconductors with broken inversion symmetry. In 
this article we focus on semiconductors with a zinc-blende 
structure where the BIA spin splitting of the electron 
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states is characterized in leading order by the Dresselhaus 

or l^gj-pn 15,21 



f^BIA(fe) 
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where 7 is a material dependent parameter. In the fol- 
lowing we will assume that the z axis is perpendicular 
to the QW plane. In asymmetric QW's SIA gives rise to 
a second contribution to spin-orbit coupling. Its leading 
term, linear in fey = (k,j;,ky,0) is frequently called the 
Rashba term.-^ It is characterized by the effective field 



f^SIA(fc||) = ^ I -krr. 

^ • 



(3) 



with a prefactor a that depends on the material parame- 
ters of the underlying semiconductor bulk material (like 
7) but in addition also on the asymmetry of the QW in 
growth direction. In first order perturbation theory for 
(001) oriented QW's the terms k^ (and powers thereof) 
in the BIA term are replaced by expectation values with 
respect to the subband function, thus one has^ 



1 
h 



kx {ky 



{kl)) 
ky {{kl) — k^) 




J^SIA(fc|| 
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Within the DP mechanism the relaxation of the com- 
ponents of S are ruled byl^ 



~ — S 7 S j 

dt ' T, ' 



(5) 



with i = z, + , — corresponding to the components of S 
along [001], [110], and [TlO], S± = l/V^{Sx ± Sy), and 
spin-relaxation ratesi^ 
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Here we have assumed an isotropic momentum scattering 
and a degenerate electron system. The rates (EJ depend 
explicitly on the carrier density n via the radius kp of the 
Fermi circle and implicitly via the self-consistent subband 
functions in the expectation values of k^ . 

In a previous evaluation of Eq. © in n-doped (001) 
grown GaAs heterostructures, Averkiev et al employed 
the triangular potential model and the variational sub- 
band function of Ref.l2^ thus assuming an infinite barrier 
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FIG. 1: Polar plot of the calculated (full line) and fitted 
(dashed line) spin splitting 2/i|r2(fe||)| in meV for fey at the 
Fermi contour in an asymmetrically doped (001) AlAs/GaAs 
heterostructure with charge density n = 3 x lO^^cm"^. 



at the interface and working in the single-band approx- 
imation. In their model, the electric field varies linearly 
with the carrier density n, which results in a oc n. Using 
Eq. lO and a tabulated value of 7 = 27eVA (appropri- 
ate for bulk GaAs) these authors calculate the relaxation 
rates for carrier densities up to 2.5 x lO^'^cm"^. (Their 
results for a smaller range of n are reproduced in Fig.[Sl) 
However, carrier densities in real AlAs/GaAs quantum 
structures with occupation of the lowest subband only 
are hmited to n < 8 x 10^^ cm^^ (see Ref. and dis- 
cussion below). Moreover, the approach of Averkiev et 
al. does not yield full self-consistency nor does it include 
nonparabolic corrections or effects due to the barrier ma- 
terial, which have to be taken into account for finite bar- 
riers. 

In order to overcome these deficiencies we perform 
fully self-consistent calculations of the subband spin split- 
ting in the multiband EFA (considering five levels or 14 
bandsi^) in application to (001) oriented AlGaAs/GaAs 
and AlGaSb/InAs QW's, n-doped on one side. Such 
calculations include in a systematic way the BIA non- 
parabolicity modified by the QW confinement and the 
SIA or Rashba term. They have successfully explained 
the anisotropic spin splitting measured with inelastic 
light-scatteringi2^ Figure ^ shows a typical polar plot 
of the calculated spin splitting in the {kx,ky) plane 
(solid contour). The striking feature is the pronounced 
anisotropy with the strongly reduced spin splitting in the 
[TlO] direct ion. From the numerically calculated spin 
splitting r2num(fe||) we extract the parameters a and 7 
by fitting r2nuni(fcF) to Eq. Q using kp = \/2nn. The 
expectation value (fc^) is evaluated using the EFA wave 
functions. The quality of this approach is demonstrated 
by the dashed contour in Fig.^ In Fig. [3 we show (k^) 
as a function of the density n together with the estimate 
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FIG. 2: Comparison of data from Ref. 13 with present results 
for the expectation value (fc^) as a function of charge density 
n for an asymmetrically doped AlAs/GaAs heterostructure. 



FIG. 4: Comparison of data from Ref. with present results 
for the BIA coefficient 7 as a function of charge density n for 
an asymmetrically doped AlAs/GaAs heterostructure. The 
dotted line indicates the value for bulk GaAs. 
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for GaA£35 and 18.5 eVA^ for AlAs. 

Recently Lau et aliSi have performed similar calcu- 
lations, yet for symmetrically doped QW's in order to 
compare with experimental data obtained from multiple 
QW's.'' Our focus here is on asymmetric single QW's and 
the interplay between SIA and BIA terms, which require 
fully self-consistent EFA calculations. 

The simpler 8-band model, which takes into account 
BIA only perturbatively, yields qualitatively the same 
spin splitting as the 14-band Hamiltonian. However, the 
magnitude of the spin splitting is about 10-20% larger in 
the 8-band model. Therefore, all results given in this 
paper are based on the more realistic 14-band model. 



III. CALCULATED SPIN RELAXATION RATES 



FIG. 3: Comparison of data from Ref. ^3 with present results 
for the Rashba coefficient q as a function of charge density 
n for an asymmetrically doped AlAs/GaAs heterostructure. 
For comparison of SIA and k linear BIA terms the dotted line 
shows the k linear part of the BIA term —'yiki). 

{kl) cx n^'^ used in Ref. Figures 01 and 0] show how 
the parameters a and 7 change as a function of n.'^'* For a 
this dependence is expected since a contains the expecta- 
tion value of the electric field. Our more realistic calcula- 
tions yield absolute values of the Rashba coefficient that 
are substantially larger than those used in Ref. |Q The 
dependence of 7 on n (which is only moderate) reflects 
higher order nonparabolicity corrections with increasing 
Fermi energy known from the bulk material:^ For our set 
of band parameters, the bulk values of 7 are 19.7 eVA'^ 



In this section we present our results for the spin- 
relaxation rates obtained from Eq. © with the param- 
eters {kl), a and 7 determined within the 14-band EFA 
calculations. This is done by assuming a momentum- 
scattering time Tp = 0.1 ps. As all rates are proportional 
to Tp they may be readily rescaled. For all calculations 
presented in this paper the range of n is chosen such that 
only the lowest (spin-split) subband is occupied. In Fig.jS] 
our results are shown for a AlAs/GaAs heterostructure 
in comparison with the data from Ref. |lj. (Note that in 
Ref. Ill To. and t_ are interchanged^^.) Both in Ref. HI 
and in our calculation, the rates and r^^ increase 
monotonously with the carrier density and are close to 
each other over the whole range. For the higher densi- 
ties (n > 0.4 X IQ-'^^cm"^) our results for rJ^ and t]^^ 
coincide with the respective data from Ref. [Tj, whereas 
our relaxation rates are generally much larger for small 
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FIG. 5: Comparison of data from Ref.^^with present results 
(marked "x") for the spin relaxation rates t^^ , , tZ^ as 
a function of charge density n for an asymmetrically doped 
AlAs/GaAs heterostructure. The range of n is chosen so that 
only the lowest (spin-split) subband is occupied. 



n. This can be attributed to the fact that in Ref. it 
was assumed that {k'^) cx n^^^, while in our calculations 
(/c^) does not tend towards zero for small n (Fig. 13) • Our 
results for are generally much smaller than those of 
Ref. ^3 and show also a pronounced nonmonotonous be- 
havior around n « 0.4 x lO-'^^cm"^, which is not present 
in Ref. ^3- A detailed explanation of the origin of this 
feature will be given below. The small oscillations below 
and above the dip at n sa 0.4 x lO^^cm^^ are numerical 
noise. 

Spin relaxation rates calculated as a function of 
the QW width L and carrier density n are shown in 
Fig- for single-side n-doped QW's with (001) orienta- 
tion for the material systems Alo.35Gao.65As/GaAs and 
Alo.3Gao.7Sb/InAs. With increasing L the asymmetric 
single-side doped QW's become more and more similar 
to a single heterostructure, and therefore the results for 
L > 200A are almost the same as for a heterostructure. 
Similar to the GaAs/AlAs system of Fig. [S] the rates 
and are close to each other in all cases, while tZ^ is 
distinctly smaller (frequently more than an order of mag- 
nitude) and shows the nonmonotonous dependence on n. 
Both the GaAs and the InAs QW's show qualitatively 
the same behavior. 

In Fig.Elwe have a dip in the density dependence of 
that is most pronounced for the InAs QW's. This feature, 
not seen in Ref. Q, appears here due to the differences 
in {k^), a and 7 (see Figs. I2I3I4|) . They cause in our 
calculations the first term in the expression defining tZ^ 
[Eq. H6b|l ] to become negative in a certain range of n and 
therefore lead to a nonmonotonous behavior. 

The origin of the dip in tZ^ can be understood more 
clearly within a picture of the DP mechanism, where we 



assume that the spin relaxation rate for the spin compo- 
nent in the direction e is only effected by those compo- 
nents of the effective field J^(fe||) that are perpendicular 
to e. In analogy with the behavior of spin relaxation rates 
for nuclear magnetic resonance^^ the spin relaxation rate 
tZ^ of the spin component in the direction [TlO] should 
be proportional to the average of the squared effective 
field components perpendicular to that direction: 

tz' (X (^{nl,„^){kp) + {nf„,,^){kp)) , (7) 

where {^'^iiQ-^)ikp) represents an average over the Fermi 
contour 

{^lio])ikF) = ^ f d^flfiio](fcf(¥')) (8a) 
^~ Jo 

and n[iiQ]{kF{ip)) denotes a projection of fl{kp{ip)) on 
the direction [110] 

fl/V2\ 

f^[iio](fei. ((/?)) = n{kF{^))- I 1/V2 j , (8b) 



cos ip 

kp{tp) = kp \ sin (/J 




(8c) 



As the projection of r2(fc||) onto the [001] direction van- 
ishes, (0^pp^j)(A:i?) is equal to zero. Hence, (ri^^^(jj)(fci?) 

is a measure of the relaxation rate r~^. In fact, by cal- 
culating this quantity as a function of n (Fig. [TJ we ob- 
tain the same behavior as for tZ^ , thus confirming the 
given interpretation. In the same line of arguments t^^ 
is essentially determined by {^'^iQ'^){kF) and, in fact, this 
quantity exhibits a monotonous dependence on n similar 
to in Figs. [SI and ini 



IV. PROPOSED MEASUREMENTS 

Our calculations indicate a large difference of the spin 
relaxation times for spins oriented in ± and z direction. 
This feature is quite universal; it appears both in GaAs- 
and InAs-based systems over a wide range of well widths 
and carrier densities. In a certain density range ("dips" 
in Fig. El the spin relaxation times can differ by more 
than two orders of magnitude. Therefore, we expect that 
this difference can be observed experimentally. 

The creation of spins oriented in z direction using opti- 
cal orientation techniques and the measurement of their 
relaxation has already been demonstrated&LSiS For the 
creation of an in-plane spin orientation a spin-injection 
technique similar to the paradigmatic spin transistor- 
can be used. Here the carrier density can be varied by 
means of a gate electrode and the spin relaxation can be 
detected via the source-drain current. 

Another way to determine in-plane spin relaxation is 
provided by measurements of the Hanle effect^S,. 
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FIG. 6: Spin relaxation rates ^, t_^^, r_ ^ as a function of charge density n for Alo.35Gao.65As/GaAs (001) (upper row) and 
Alo.3Gao.7Sb/InAs (001) (lower row) quantum wells with widths L — 50, 100, 200 A and for a simple heterostructure. The 
range of n is chosen so that only the lowest (spin split) subband is occupied. 
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FIG. 7: Mean squared transverse field (fi^j^j^oi)^^^) ^'^ ^ func- 
tion of n for the AlGaAs/GaAs systems of Fig. |S| 



We propose a second scheme for the creation of in- 
plane excess spins and the measurement of their relax- 
ation time: In a first step, spins oriented in z direction 
are created with conventional optical orientation tech- 
niques. These can be rotated into the plane of the QW 
by means of a so-called tipping pulse, as demonstrated by 
Gupta et alSL After a certain delay, the spins are rotated 
back into the z direction, so that the magnitude of the 
remaining spin polarization can be measured. By chang- 
ing the delay time the spin relaxation for an arbitrary 
in-plane direction can be analyzed. A substantial differ- 
ence between the spin relaxation times can be expected 
for spins oriented along the in-plane directions [110] and 



We propose that such a system can be used as a "spin 
storage" : The rotation of spins from the z direction into 
the [110] direction effectively stores the spins because the 
relaxation time for the latter direction is larger by orders 
of magnitude. 
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V. CONCLUSION 

Starting from a realistic model for the anisotropic 
spin splitting in GaAs/AlGaAs and InAs/AlGaSb quan- 
tum structures we calculate the parameters which deter- 
mine the longitudinal (r~^) and transverse (r^^) spin- 
relaxation rates and determine their dependence on QW 
width and carrier concentration. We quantify the previ- 
ously predictedii giant anisotropy of spin relaxation by 
our more realistic calculation of the spin splitting. In 
contrast with previous investigations our accurate calcu- 
lations predict a pronounced nonmonotonous behavior of 
rZ^in), which can be understood by ascribing the spin 
relaxation to the transverse components of the effective 
magnetic field. Future work will include calculations for 



different material systems, growth directions and doping 
profiles. For comparison with room-temperature exper- 
iments the temperature dependence and the effects of 
different scattering mechanisms have to be considered. 
Furthermore, the microscopic interface asymmetry, that 
has been identified as an additional mechanism of spin 
splitting, "^^ will be taken into account. 
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